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ABSTRACT 

The brightness of gamma-ray burst (GRB) afterglows and their occurrence in young, blue galaxies make 
them excellent probes to study star forming regions in the distant Universe. We here elucidate dust extinction 
properties in the early Universe through the analysis of the afterglows of all known z > 6 GRBs: GRB 090423, 
080913 and 050904, at z = 8.2, 6.69, and 6.295, respectively. We gather all available optical and near-infrared 
photometry, spectroscopy and X-ray data to construct spectral energy distributions (SEDs) at multiple epochs. 
We then fit the SEDs at all epochs with a dust-attenuated power-law or broken power-law. We find no evidence 
for dust extinction in GRB 050904 and GRB 090423, with possible evidence for a low level of extinction in 
GRB 080913. We compare the high redshift GRBs to a sample of lower redshift GRB extinctions and find a 
lack of even moderately extinguished events [Ay ~ 0.3) above z > 4. In spite of the biased selection and small 
number statistics, this result hints at a decrease in dust content in star-forming environments at high redshifts. 
Subject headings: dark ages, reionization, first stars - dust, extinction - galaxies: high-redshift - gamma-ray 
burst: individual (GRB 090423, 080913, and 050904) 



1. INTRODUCTION 

Dust formation in the early universe is a hotly debated topic. 
The properties and quantity of dust at these times are es- 
sential not only to our understanding of how dust forms and 
evolves, but also to models of star-formation and the appear- 
ance of the first generations of stars and galaxies. Currently 
we extrapolate dust properties in the local group to very dif- 
ferent environments that existed less than a billion years af- 
ter the Big Bang. Long-duration gamma-ray bursts (GRBs) , 
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associated with the deaths of massive stars (IWoosle 

Galama et al. 1998; iHiorth et al.1 [2 003: St anek et al. 

Malesani et al.i i2004t ICampana et al.i i2006) . have extremely 
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brigh t and spectrally simple afterglows (Fy oc v"^) jSari et alj 
119981: iGranot & Saril 120021) . and so provide an opportunity 
to obtain not only effective reddening, but also absolute ex- 
tinctions from distant galaxies (e.g . IWatson et al.l l2006al : 
iSchadv et al.ll20ldtlZafar et al.ll201 lb . The recent tremendous 
advances in thi s subject were enab led by the dedicated GRB 
satellite. Swift jGehrels et al.ll2004l) . with its precise localiza- 
tion and fast responding X-ray and UV/optical telescopes, de- 
tecting the afterglow within only a few tens of seconds after 
the y-ray trigger The data from Swift have provided a very 
large sample of afterglows virtually complete in X-ray detec- 
tions, and with a high completeness in the optical. It has also 
for the first time, enabled the discovery of GRBs at redshifts 
greater than or comparable to the most distant know n galax- 
ies and quasars (Q SO) (see Kawai et al. ,2006 : .Greiner et alj 
l2009tlTanvir et alj|20 09: Salvaterra et al.1 l2009l) . 

In this paper we collect for the first time all the available 
afterglow photometry and spectra from near-infrared (NIR) to 
X-ray at different epochs for all GRBs at z > 6: GRB 050904 
(z = 6.295), GRB 080913 (z = 6.69) and GRB 090423 
(z ~ 8.2). Since observations of these high-z GRBs are be- 
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ing made in the rest frame UV, they are in principle relatively 
sensitive to even small amounts of dust. We extract the SEDs 
at multiple epochs for all three GRBs and jointly fit the X-ray 
to NIR data to determine the properties of dust at z > 6. In §2 
we describe multi-wavelength observations of the afterglows 
at different epochs carried out with different instruments. In 
§3 we present our results from the SED fitting. Based on 
our results, in §4 we discuss the origin of dust at high red- 
shift and compare to the first large spectroscopic sample of 
GRB-derived extinction curves, the first spectroscopic sam- 
ple of absolute extinction curves outside the local group. In 
§5 we provide our conclusions. 

For all parameters we quote uncertainties at the 68% confi- 
dence level. We derive 3cr upper limits for the cases where 
detection is less than 2cr significant. A cosmology where 
Hq = 72kms"' Mpc"', Qa - 0.73 andfim = 0.27 is assumed 
throughout. 

2. MULTI-WAVELENGTH DATA 

The X-ray to NIR SED o f GRB 050904 wa s analysed in 
detail in a previous work (see lZafar et al.ll20Tol) . 

On 2008 September 13 at 06:48:33.6 UT Swift-XRT began 
observing GRB 080913. Telescopes at the European Southern 
Observatory (ESO), Gemini Observatories, and the Subaru all 
observ ed the field of GRB 080913 in different optical and NIR 
bands (iGreiner et al.ll2009.) . On 2008 September 13, an opti- 
cal spectrum of the afterglow of GRB 080913 was secured at 
the Very Large Telescope (VLT) using the FOcal Reducer and 
low dispersion Sp e ctrograph 2 (FOR S2) (IGreiner et al.ll2009l : 
iFvnbo et al.1 120091: iPatel et al.ll20Tol) . The spectrum shows a 
sharp break at ~ 9400 A supporting the detection of Lya ab- 
sorption at z = 6.69. We use this spectrum in addition to pho- 
tometric dat a in the J, H and K bands from the lightcurves 
presented in IGreiner et al.l (l2009h . We corrected the data for 
extinction in the Milky Way (MW) of E(B - V) = 0.043 mag 
along the line of sight to this burst (Schlegel et al. 1998). 

At 07:56:32 UT on 2009 April 23, Swift-XKT began obser- 
vations of GRB 090423. The afterglow of GRB 090423 was 
observed in optical and NIR bands, wit h non-detection in g riz' 
and y-band due to the high redshift dTanvir et al.l 120091) . A 
NIR spectrum was obtained with the VLT Infrared Spectrom- 
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Table 1 

Best fit parameters of the afterglow SEDs at different epochs with Small 
Magellanic Cloud (SMC) type extinction. The provided upper limits are 
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GRB 


Epoch 


Av 
(mag) 


/Si 


P2 


log Vtreak 

Hz 


050904 


0.47d 
1.25d 
3.4d 


< 0.17 

< 0.08 
<0.16 


1.23 ±0.08 
1.17 ±0.51 

1.24 ±0.07 






080913 


10.5m 
26.9m 
1.88hr 
4.78d 


0.14 ±0.06 
0.22 ±0.11 
0.12 ±0.03 
< 0.40 


0.78 ± 0.02 
0.65 ± 0.04 
0.79 ± 0.03 
0.94 ± 0.04 






090423 


16.7hr 
1.65d 


< 0.10 

< 0.09 


0.45 ±0.11 
0.68 ±0.13 


0.95 ±0.11 
1.18±0.13 


15.1 ±0.7 
16.5 ± 0.9 



eter And Array Camera (ISAAC) starting about 17.5 hours 
after the burst trigger A redshift of ~ 8.2 was estimated 
from a break observed around 1.14/im due to Lyo- absorp- 
tion by the neutral intergalactic medium (IGM). Another spec- 
trum of the afterglow was obtained 40 hours after the burst 
with the VLT's Spectrograph using INtegral Field Observa- 
tions in the Near Infrared (Sl NFONl) and confi rms the anal- 
ysis of the ISAAC spectrum jTanvir et al.ll2009l) . The finding 
was also confirmed with a spectrum obtained with Telesco- 
pio Nazionale Galileo (TNG) at ~ 14 hours after the burst 
( iSalvaterra et alj|2009l). We use the /, H and K band pho- 
tometrv of iTanvir et alj (|2009). The data wer e corrected for 
Galac tic extinction of E{B -V) - 0.029 mag dSchlegel et alj 
[19981) . 

For all three GRB afterglows the Swift XRT data were avail- 
able and downloaded from the Swift data archive and reduced 
using HEAsoft (version 6.10). Spectra and lightcurves were 
extracted in the 0.3-10.0 keV energy range. The X-ray spec- 
tra at different epochs were analyzed using the latest associ- 
ated calibration files. All X-ray spectra were obtained near 
the time of optic al/NlR spectra or relevant available photom- 
etry. To get an approximate X-ray flux level at each epoch, 
the afterglow lightcurves were obtained from the Swift XRT 
GRB light curve repository at th e UK Swift Science Data 
Centr^ created as described in lEvans et al] ( 120071 120101) 
and fitted by assuming a smoothly decay broken power-law 
(iBeuermann et al.lll999h . Considering the photon weighted 
mean time, the X-ray spectra were then normalized to the rel- 
evant SED time by using the lightcurve fit. 

The XMM-Newton observations of the afterglows of 
GRB 080913 and GRB 090423 were performed at ~ 4.4 and 
~ 2.2 days respectively. The data were collected from the 
XMM.-Newton science archive (XSAfl The data were re- 
duced in 0.3-10.0 keV band using Science Analysis System 
(SAS, version 10.0.0). The spectrum of GRB 080913 was 
fitted using a single power-law with a best fit photon index 
of r = 2.4 ± 0.4 and frozen for Galactic X- ray absorption 
of 3.2 X 10^0 cm - (using the nH FTOOL; iKalberla et af] 
12005b . The best fit photon index for GRB 090423 is F = 
2.4 + 0.3 with the fixe d Galactic absorption of 2.9 x 10^" cm - 
(IKalberla et aUlMBb . 

Parts of spectra with a restframe wavelength shortward of 
Lya (i.e. A < 1216A) are not included in the SED analysis to 

' www , swift . ac . uk/xrt_curves| 
^ http: //xmm. esac .esa. int/xsa7| 
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Figure 1. NIR to X-ray SED of the afterglow of GRB 080913 at 10.5 
minutes (open red triangles), 26.85 minutes (blue circles), 1.88 hours 
(green squares) and 4.78 days (black triangles). Yellow stars correspond to 
^yBA-Newton data at 4.4 days after the burst. The grey curve is the optical 
spectrum of the afterglow taken at fo + 1.88 hours. The solid lines represent 
the best-fit to the data; dashed lines illustrate the unextinguished power-law 
from the best fit. Being close in time, the X-ray specfi'a for the first two 
epochs are superimposed. 

avoid absorption caused by the neutral IGM. The absorption 
lines arising from ionic species were removed from the spec- 
trum of GRB 050904. Uncertain ty in the Galactic e xtinction 
does not affect our re sults. The iDutra et alj (120031) redden- 
ing maps confirm the ISchlegel et al.l (119981) Galactic extinc- 
tion maps up t o E{B — V) - 0.25 mag. Assuming an uncer- 
tainty of 15% (ISchlafly et al.ll2010l) even for the case with the 
largest Galactic extinction (z'-band, GRB 050904), this cor- 
responds to an uncertainty in the Galactic reddening value of 
< 0.015 mag. This indicates that the uncertainty due to the 
Galactic extinction correction in the NIR is always smaller 
than our statistical uncertainties. 

3. RESULTS 

SEDs of the three GRB afterglows were obtained at multi- 
ple epochs. The data were fitted with dust-attenuated power- 
laws or broken power-laws. In addition to our standard ex- 
tincti on model (SMC type, see lSchadv et al.l2010tlZafar et alj 
1201 ll) . w e also fitted the data with the extinction curve in- 
ferred by iMaiolino et alJ (|20()4|) to explain the spectrum of a 
QSO at z = 6.2 (although see iGallerani et al ...20 1 tt) . 

3.1. Multi epoch SED of GRB 050904 

The SED of this burst was examined in detail in lZafar et alj 
(I20T0I) . where the SED of the afterglow of GRB 050904 is 
constructed at 0.47, 1.25 and 3.4 days after the burst. We 
use the results from the analysis here. The results of the 
SMC-extinction fits are reported in Table. [T] for complete- 
ness. The fits using supernova (SN)-origin extinction from 
IMaiolino et alj ( 120041) suggest that there is no evidence in this 
burst for this type of extinction. Th e SED of GRB 050904 ha s 
been fitted by HaisIicsLaLl (120061) : iTaghaferri et"an (120051) : 
iKann et aP (l2010i) finding no dust extinction. 

3.2. Multi epoch SED of GRB 080913 

The SED of the afterglow of GRB 080913 is extracted 
at four epochs i.e. 10.5 and 26.9 minutes, 1.88 hours and 
4.78 days after the burst (see Fig. [Til. To get precise photome- 
try in the z-band which is affected by strong Lya absorption. 
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GRB 080913 (t=1.88hr) 
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Figure!. NIR and optical SED of the afterglow of GRB 080913 at 1.88 
hours after the trigger. The grey curve represents the optical spectrum. 
The solid black curve corresponds to the best fit SMC-dust-attenuated 
power-law to the full NIR-X-ray SED, while the red line shows the best 
fit power-law to the NIR-optical data only with no extinction. Inset: Zoom 
of the optical spectrum and fits. The spectrum has been median filtered for 
display purposes. 



we performed spectro-photometric analysis at 1.88 hours us- 
ing the effective transmission function of the FORS2 z-Gunn 
filter. We constructed SEDs of the st ars in the field using 
the Sloan Digital Sky Survey (SDSS; iFukudta et allll996b 
and Two Micron All Sky Survey (2MASS: ISkrutskie et all 
|2006) catalogues. We derived the count-to-flux conversion 
factor and estimated the afterglow band-integrated flux from 
its measured counts. We computed the flux density at A^est - 
1250 A using the procedure described in IZafar et al.l (120101) . 
At 1.88 hour the flux density at 1250 A was 2.66 + 0.63 /iJy. 
The NIR to X-ray SED at all epochs is well fitted with a sin- 
gle power-law and SMC type extinction {x~/dof - 47/81). 
The combined fit for dust extinction at all epochs resulted in 
marginal dust reddening of Ay = 0.12 + 0.03. Fitting each 
epoch individually, we find dust extinction significant at the 
2cr level at 10.5 and 26.9 min, at 4cr at 1.88 hours and Icr 
at 4.78 days (see Tab le [T]i. We also fitted the data using the 
iMaiolino et al.l (l2004h extinction curve, resulting in a signif- 
icantly worse fit (x^/dof = 69/81) and zero extinction. We 
also fitted the SED with Milky Way and Large Magellanic 
Cloud extinction curves but found no statistically significant 
improvement in the fit compared to the SMC, and we found 
no evidence for a 2175 A bump. As the simpler model, we 
used the SMC to parameterize the extinction here and to pa- 
rameterize the upper limit in GRB 090423 below. 

It is worth noting that the X-ray lightcurve for the afterglow 
of GRB 080913 shows some evidence of flaring, possibly re- 
sultin g in harder spectra (see IE vans et al.ll2007l; iGreiner et al.l 
l2009b . If the SED of the afterglow of GRB 080913 is domi- 
nated by two separate components in optical-NIR and X-ray 
wavelengths, as observed at early times for GRB 050904, we 
should fit the optical-NIR data independent of the X-rays. Us- 
ing this method, we obtain results consistent with a single 
power-law and no dust extinction, with a spectral slope of 
/3 - 1.2 + 0.1 and an upper limit of Ay < 09 (se e Fig. |2]i. 
Previously iGreineret al. (2009); K ann et al.1 (120101) fitted the 
optical-NIR SED of GRB 080913 and found no dust extinc- 
tion. 
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Figure 3. The NIR to X-ray SED of the afterglow of GRB 090423 at 16.7 
hours (black triangles) and 1.645 days (blue squares). The corresponding 
fines show the best fits to the data. The ISAAC and SINFONI spectra are 
shown with the grey and cyan curves respectively. Red stars represent the 
'KMM-Newton data taken at 2.21 days after the trigger. 



3.3. Multi epoch SED of GRB 090423 

The SED of the afterglow of GRB 090423 was constructed 
at 16.7 hours and 1.654 days (see Fig.[3l). The VLT/ISAAC 
and SINFONI spectra, NIR photometry and X-ray spectra 
were fitted at both epochs. The data from both epochs were 
well fit with no dust using a broken power-law. The estimated 
3<T upper limit on the Ay with SMC extinction is 0.10 and 
0.09 mag at 16.7 hours and 1.654 epochs respectively. 

4. DISCUSSION 

Two apparently conflicting themes appear in studies of dust 
at high redshift. The first is the expected decline in extinc- 
tion as we get to early in the Universe, where galaxies are 
younger and dust has less time to form. In the local Uni- 
verse low and intermediate mass asymptotic giant branch 
(AGB) stars are believed to be the principal sources of in- 
terstellar dust, requiring > 1 Gyr to ev olve off the main se- 
quence to produc e dust (Gehrz 1989; Mor gan & Edmundsl 
120031: lMarchenkdl2006; Dweketal. 2007). Because of the 
time required to evolve to the giant branch, it has been posited 
that AGB stars should n ot contribute si gnificantly to dust a t 
z > 6 dDwek et al.l2007l). Spectrosco pic dDouglas et al.l2010l) 
and photometric (IVerma et al.ll2007l) studies of z ~ 5 lyman- 
break galaxies seem to confirm the general trend of little dust 
at high redshifts, showing low extinctions, and the sample 
of z-band dropout ga laxies (z ~ 7) shows little or no dust 
dGonzalez et al.ll2010l) . 

However, the second, at first paradoxical, theme is the 
very large dust masses discovered in quasar h ost galaxies at 
z ~ 6 from mm and sub-mm observations (e.g. Bertoldi et alj 



2003l;ICarilli et alj,200kiBeelen et a 


.,2006 


^,Wang et al..200l 


Michalowski et al.ll2010l; iGaU et al.1 


2011b 


), and the mystery 



surrounding the origin of such large quantities of dust within 
a billion years of the Big Bang. It has been suggested that dust 
produced on short time scales could be produced in the ejecta 
of Type II SNe (Kozasaetal. 1991; Todini & Ferrara 20011; 
Morgan & Edmunds 2003; Nozawaetal. 2 0(j3l: lDwek et alJ 
2007i; iMarchenka,2006; .Hirashita et al.ii2005ir although ob- 



servations of local SNe do not suppo r t this hypothes i s (e.^ 
iMeikle et al.l [20071; iSmith et al.ll2008l) . iValiante et al. l (ilOOl 
estimated that under certain conditions AGB stars could pro- 
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duce dust at earlier epochs, beginning to dominate b etween 
150 and 500 Myr However iMichalowski et al.l ( 120101) found 
that given the stellar masses in high redshift QSO hosts, AGB 
stars could not be responsible for the l arge masses of dust 
observed i n these objects at z > 5 (see Che rchneff & Dwe3 
[2009; Pi pino et al.l 1201 It iGaU et al. 2011a). and su ggested 
rapid dust growth in the ISM as an alternative (see also lDrain^ 
120091) . The modeling of Lyman alpha emitters at z ~ 6 sug- 
gests that the largest galaxies are the most dust enri ched and 
that dust content decreases with increasing redshift (iLai et alJ 
12007; Finkelstein et al. 2009; Dayal et al.'2010). In spite of 
the emphasis in the literature on these QSO hosts, it may 
be that such very massive systems, being exceptional and 
rare, tell us little about typical star forming environments at 
very high redshift. And while, young star-forming galax- 
ies at high redshift with strong nebular emission appear to 
be somewhat du st-obscured in spite of their blu e continuua 
(ISchaerer (feleBarros 20Tot IWatson et al.ll2010l) . the quan- 
tities of dust are not large, consistent with the idea that in 
the early universe, most of the star-formation is likely to 
be less obscured than at later times. GRBs, associated as 
they are with the deaths of massive stars, usually occur in 
young, blue and sub-lumino us galaxies with a high specific 
star formation ra te (Le Floc'h et al.1l2003 ; Courtv et al1l2004l: 
Christensen et alji2004; f rochaska et al.ll2004; .Savaglio et alj 
2009HCastro Ceron et alJl2oTol) ' and are therefore likely to be 
good probes of star-forming regions, especially at high red- 
shift. 

In this work, we have determined that all the identified 
z > 6 GRBs have low or negligible extinctions. By compar- 
ing these r esults with t he firs t spectroscopic GRB extinction 
sample of IZafar et al.l (120111) . we find that at high-redshifts 
(z > 4), GRB sightlines appear significantly less extinguished 
(see Fig. |4]i: the majority of reddened low-redshift GRBs are 
observed with Ay ~ 0.3 mag. While the numbers are small, 
the absence of any GRBs at z > 4 with extinctions in this 
range is striking. The first explanation to examine is observa- 
tional bias, where higher-redshift GRBs will be observed far- 
ther into the restframe UV where the extinction is likely to be 
more severe, possibly preventing spectra from being obtained 
or the afterglow from being detected. For the z > 6 GRBs, a 
restframe Ay = 0.3 corresponds to 1-1.5 magnitudes of ex- 
tinction in the observed /-band and 0.5-0.8 magnitudes in the 
observed K, for an SMC extincti on curve. Considering th e 
p hotometric observa ti ons given inlRuiz-Velasco et al.l (120071); 

faislip et al.' (2006); Ta diaferri et al.l (120051) : iGreiner et all 
009); TanyireLaL(2009t), we found that even with Ay = 0.3 
the GRBs at z > 6 could have been detected with rela- 
tive ease and a photometric redshift obtained. In the case 
of GRB 050904, it seems likely that a spectroscopic redshift 
would have been obtained even with restframe Ay - 0.3. 
This work therefore hints at a low extinction of most z > 4 
bursts, since it seems that our z > 6 bursts would have been 
found even with Ay = 0.3. Therefore we have looked at all 
known z > 5 GRB afterglows in the literature. There are four 
bursts with firm estimates of z > 5 known. Three of them 
have estimates of their d ust extinction availabl e. The four 
bursts are: GRB 050814 ( Jakobss on eta^l2006^ at z ~ 5.3, 
GRB 060522 (Cenkoetal. 2006) at z = 5.11, GRB 060927 
dRuiz-Velasco et al. 2007; Zafar et an i2011l) at z = 5.46 and 
GRB 071025 (.Perlev et al...201Q) at z ~ 5. iJakobsson etal] 
(|2006) obtained a photometric redshift of the afterglow of 
GRB 050814 and the SED is consistent wi_th no dust redden- 
ing (P. Jakobsson, private comm. 2010). ( iCenko et al...2006h 
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Figure 4. Extinction in GRB afterglow versus redshift. The three high 
redshift GRBs (z > 6) are fro m this work wh ile the 42 low redshift GRBs are 
from a sample studied in Z afar et al J 1201 If) . There seems to be a dearth of 
GRB afterglows at z > 4 with Ay ~ 0.3. 



give a spectroscopic redshift for GRB 060522 based on the 
Lya break, however no e stimate of the extin ction is available 
so far in the literature. iRuiz-Velasco et all (12007) obtained 
the spectroscopic redshift of the afterglow of GRB 060927. 
We have studied the SE D of the afterglow of this burst in our 
spectroscopic sample of lZafar et al.l (1201 II) and it is consistent 
with Ay - Q (see Fig. |4|. Moderate extinction is claimed in 
the case of GRB 071025 (iPerley et al. 2010) , with extinction 
corresponding to Ay ~ 0.5 (conv erting A^nnn k to Ay u sing 
the SMC extinction curve ratio). A lMaioUno et alj ( 120041) ex- 
tinction curve is required to fit the SED. With this extinction, 
GRB 071025 has a high enough extinction to begin to fill in 
the gap at high redshift around Ay ~ 0.3, though it is perhaps 
notable that it has the lowest redshift of all the bursts exam- 
ined here. It is more anomalous in the sense that it has an ap- 
parently unique extinction curve since reanalyses have shown 
that no other QSO or GRB any long er require such a p eculiar 
extinction (see lGallerame t al. 2008: lZafar et al.ll20ld) . 

The issues of bias and low number statistics are clearly im- 
portant and we cannot resolve them with this sample. The 
fact that we infer very little dust in the afterglows studied here 
does not imply that dusty environments do not exist along 
sightlines to z > 6 GRBs, but it suggests that there is less dust 
in highly star-forming regions at z > 4. This result is also po- 
tentially interesting in the context of dark bursts - suggesting 
that highly-extinguished bursts may not exist at high redshift, 
and that bright, very high redshift GRBs will be detectable as 
long as they are observed in the NIR. 

Finally, the comparison to the column densities of met- 
als is striking. In the cases of both GRB 050904 and 
GRB 090423, column densities of Nh,x ~ 3 x lO^^ ^nd 
1 X lO^-' are claimed (assuming solar abundances and neu- 
tral elements, |W atson et al. 2006R ICampana et al.l 1201 ll : 
iGendre et al.ll2006i: iTanvir et al.ii20()9l : Starling et al. in prep. 
201 1) from the soft X-ray absorption. The dust-to-metals ra- 
tios in both cases are at most a few percent of the value in the 
local gr oup. Given th at the metal s are likely to be highly- 
ionized dSchadv etalj|20i 0; Wats on et all 120071) . the metal 
column must be substantially higher than found using neu- 
tral medium models. This could drive the dust-to-metals ratio 
down below a percent of local group values. While it is estab- 
lished that GRB environments can have low dust-to-metals ra- 
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tips (lEUasdottir et al.l2009l: I Watson et al.l2006aHSchadv et alj 
l20icil) . these high-redshift objects are extreme. Such extreme 
ratios could be related to the more effective dest ruction of dust 
due to the high luminosity of these bursts (e.g. lFruchter et all 
120011) : the effects of dust destruction may be discernible not 
only in very early colour changes in the GRB prompt phase 
UV emission, but also in the extinction curve at later times 
jPerna & Lazzatil [20021) . A small dust-to-metals ratio could 
however, also be due to something more fundamental, such 
as the rapid production of metals in SNe but delayed forma- 
tion and growth of the dust. For the metals to form before 
the dust suggests at least several million years between the 
formation of metals and the onset of significant dust growth. 
Obs ervations of a larg er sample may help to clarify this issue 
(see lZafaretai]|20m . 

5. CONCLUSIONS 

We have investigated cosmic dust at z > 6 using the after- 
glows of the highest redshift GRBs known. We performed 
multi-epoch NIR-X-ray SED analysis of GRB 050904, 
080913, and 090423. We infer from our analysis that there 
is no evidence of dust in GRB 050904 and GRB 090423. We 
find possible evidence of low extinction in GRB 080913 con- 
sistent with an SMC extinction curve, however the optical 
spectrum and NIR photometry alone are consistent with a 
power-law with no dust extinction. In no case did we find 
evidence for any oth er extinction curve (e.g. the 'SN-origin' 
curve of iMaiolino et al. 2004). Comparing to a much larger 
spectroscopic sample of GRB extinctions, we find a distinct 
absence of extinction at high redshifts. While the high red- 
shift sample is very small and some of the effect may be ex- 
plained by restframe UV-selection bias, this is not the whole 
picture, as all the z > 6 GRBs would have been detected at 
the typical Ay of lower redshift reddened bursts. This hints 
that there is less dust along sightlines to the highest redshift 
GRBs, indicating less dust in the early Universe. From these 
results, we also infer an extremely low dust-to-metals ratio in 
GRBs at high redshift, suggestive either of efficient dust de- 
struction, or a delay of at least several million years between 
the formation of metals and the formation and growth of dust. 
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